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Abstract. In this paper, we obtain a blow up criterion for strong solutions to the 
3-D compressible Naveri-Stokes equations just in terms of the gradient of the velocity, 
similar to the Beal-Kato-Majda criterion for the ideal incompressible flow. The key 
ingredients in our analysis are the a priori super-norm estimate of the momentum by 
a Moser-iteration and an estimate of the space-time square mean of the gradient of 
the density. In addition, initial vacuum is allowed in our case. 

1. INTRODUCTION 

We shall study the following isentropic compressible Navier — Stokes equations in 
3 — V case: 

'dp T , . 
— + div(/)u) = 0, 

+ div(pu u) - ^Au - + A)V(divu) + VP(p) = 

Where p, u, P denotes the density, velocity and pressure respectively. The pressure- 
density state equation is given by 

P{p) = ap^ (a>0,7>l) (1.2) 

p and A are shear viscosity and bulk viscosity respectively satisfying the physical con- 
dition: 

p>Q,\ + ^p>Q (1.3) 

Lions [1] [2], Feireisl [3][TT] et. established the global existence of weak solutions to 
the problem (1.1) — (1-3), where vacuum is allowed initially. The global existence to the 
compressible Navier-Stokes equations is obtained by Matsumura|18jand Nishida under 
the condition that the initial data is a small perturbation of a non-vacuum constant. 
It is also shown by Xin^l] that there is no global in time regular solutions in E? 
to the compressible Naiver-Stokes equations provided the initial density is compactly 
supported. 
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There are many results concerning the existence of strong solutions to the Navier- 
Stokes equations, only local existence results have been established, see [E], [E] , [TT] , [20] . 
V.A.Solonnikov proved in [TH] that for pressure laws and initial data satisfies for 
some q > N, 

0<m< po{x) <M <oo, and po £ W^^'^(r^) (1.4) 

there exists a local unique strong solution (p, u) to (1.5) — (1.6) for periodic data, such 
that 

p G L°°(0, T; W^'^{T^)), pt G L'^((0, T) x T^) 

(1.6) 

u G L'?(0, T; W^'''{T^)), ut £ L^O, T) x T^f 
Later, it was shown in [T5| that if is either a bounded domain or the whole space, 
the initial data po and uq satisfy 

< Po G VF^'^(O), UQ G i7o^(0) n (1.7) 

for some q G (3, oo) and the compatibility condition: 

- pAtio - (A + p)Vdivuo + VP(po) = Po^^g for some g G l2(0) (1.8) 

then there exists a positive time Ti G (0, oo) and a unique strong solution (p, u) to the 
isentropic problem, such that 

p G C([0,ri];t^i'*(17)), Pt G C([0,ri];L^«(J7)) 

u G C([0,Ti];D^ nL»2(^)) ^2.2(0,ri;D2,<?o(5^)) (i.g) 

Ut G l2(0, Ti; i?i(^^)) G L°°(0, Ti; l2(17)) 

Furthermore, one has the following blow-up criterion: if T* is the maximal time of 
existence of the strong solution (p, u) and T* < oo, then 

sup (||p|liyi,<7o + \\u\\d}) = oo (1.10) 

where q^ = min(6, g). 

Throughout this paper, we use the following notations for the standard homogeneous 
and inhomogeneous Sobolev spaces. 

D^^-'in) = {u£ Ll^in) : llV'^uWir < oo}. 
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Dq = {u G L^(r2) : ||Vii||^2 < Qo and u = on dCl}, 

Hl=L^r\Dl, \\u\\Dk.r = \\V''u\\Lr 

Very recently, Fan and Jiang [H] proved a blow-up criterion for such strong solutions, 
i.e, when 7/i > 9A, 

lim { sup IIpIIloo + / (||p||tyi.,o + ||V/9||^2)di} = oo (1.11) 

T^T* o<T<T* Jo 

Here they only require a sufficient regularity of density p to admit the global existence 
of strong solutions, as (1-11) revealed. 
In this paper, we assume that 

^ + A > 0, = 2, = 

(1.12) 

^ + A = 0, N = 3, ncR^ 
Here and thereafter C always denotes a generic constant depending only on fi, T and 
initial data. 

For the initial boundary value problem, we have the following result: 

Theorem 1.1. Let Q be either a three dimensional bounded domain or two dimensional 
torus. Qt = (0,T) X 0. Assume that the initial data satisfy (1.7) — (1.8) and (1.12) 
holds. Let (p, u) be a strong solution of the problem (1.1) — (1-3) satisfying the regularity 
(1.9). If T* < oo is the maximal time of existence, then 

lim / ||Vu||^oo(Q)(it = oo (1-13) 

J 

In the case of initial value problem, it holds that 

Theorem 1.2. Let = B? . Assume that the initial data satisfy 

po e H\R^) n W^^^iR^), uo G Dl{R^) n D\R^) (1.14) 

for some q with 3 < q < oo and the compatibility condition (1.8). Let {p,u) be a strong 
solution to the problem (1.1) — (1-3), and satisfy 

p e C{[0,Ti],H^{R^) nW^''i''{R^)), pt e C{[0,Ti],L^{R^)nL'"'{R^)) 

u G c([o,ri],z)^(i?3) nL>2(ii3)) nL2(o,ri;L>2''?"(i?3)) (i.i5) 

ut G L^{Q,Ti;Dl), ^put G L°°(0,ri;L2(i?3)) 
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where qo = min(6,g). If T* < oo is the maximal time of existence, then 




■T 



lim 



(1.16) 



Remark 1.1 The blow up criterion (1.10) involves both the density and velocity. It 
may be natural to expect the higher regularity of velocity if density is regular enough. 
(1.11) shows that sufficient regularity of the gradient of density indeed guarantees the 
global existence of strong solutions. The main difficulty in our case is to control the 
gradient of density, which is not a priori known and coupled with the second derivative 
of velocity. 

We develop some new estimates under the condition (2.1). In fact, the key estimates 
in our analysis are both bound of pu and L°°(0, T; L^(r2)) norm of Vp. The Super- 
norm estimate for the momentum is obtained by a Moser-iteration based on the a priori 
energy bounds motivated by an analysis in [T^]. To control the L°°(0, T; L^(r2)) norm 
of V/3, our key observation is that the space-time square mean of the convection term 
F = put + pu ■ Vn is controlled by that of V/> (see Lemma 2.3). This, in turn, gives the 
desired L~(0, T; ^^(J])) estimate on Vp, and thus the L'^{0,T; H'^{n)) of u. Then the 
higher order regularity can be obtained by using the equations and the compatibility 
condition (1.8). 

Remark 1.2 There are many results concerning blow-up criteria of the incompress- 
ible flow. In their well-known paper 0], Beal-Kato-Majda established a blow-up cri- 
terion for the incompressible Euler equations. One can get global smooth solution if 
\\uj\\Loadt is bounded. It's worth noting that only the vorticity uj plays an impor- 
tant role in the global existence of smooth solutions. Moreover, as pointed out by 
Constantin[9], the solution is smooth if and only if /q^ ||((Vti).^) • is bounded, 

where ^ is the unit vector in the direction of uj. It turns out that the solution be- 
comes smooth either the asymmetric or symmetric part of Vn is controlled. Later, 
Constantin[7], Fefferman and Majda showed a sufficient geometric condition to control 
the breakdown of smooth solutions of incompressible Euler involving the Lipschitz reg- 
ularity of the direction of the vorticity. It is also shown by Constantin[8] and Fefferman 
that the solution of incompressible Navier-Stokes equations is smooth if the direction 
of vorticity is well behaved. 
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Recently, in[5], assuming that the added stress tensor is given in a proper form, and 
using an idea of J.-Y. Chemin and N. Masmoudi [6], Constantin, P. and Fefferman, 
C, Titi, E. S. and Zarnescu, A obtain a logarithmic bound for \\\7u\\Loodt. to con- 
clude that the solution to Navier-Stokes-Fokker-Planck system exists for all time and 
is smooth. 

In our paper, we establish a similar criterion to Beal-Kato-Majda. Our blow up 
criteria involve both the symmetric and asymmetric part of Vn, as compressibility and 
vorticity are two key issues in the formation of singularity of compressible Navier-Stokes. 

Remark 1.3 The condition (1.12) is assumed to obtain L°°(0, T; L°°(r2)) norm of 
the momentum pu by a Moser-Iteration. However, in 2 — D periodic case, using an 
estimate motivated by Desjardin[TO], we can show (1.13) holds for natural physical 
constraint: ^ + A > 0. 

Remark 1.4 We will study the blow up criteria for smooth solution of compressible 
Navier-Stokes in another paper [13]. 

2. Proof of Theorem 1.1 

Let {p,u) be a strong solution to the problem (1.1) — (1-3). We assume that the 
opposite holds, i.e 

_^lirn^ J ||Vn||ioo(Q)(it < C < oo (2.1) 
First, the standard energy estimate yields 

sup ||p^/\(t) 11^2 + / \\u\\jj^dt<C, 0<T<T* (2.2) 

0<t<T Jo 

By assumption (2.1) and the conservation of mass, the bounds of density follows 
immediately, 

Lemma 2.1. Assume that 

[ \\divu\\Loodt<C, 0<T<T* (2.3) 

then 

||pIIl-{q^) < o<r<r* (2.4) 

Proof. It follows from the conservation of mass that for Vf? > 1, 

dtip") + div(p^u) + (q - l)p^divu = (2.5) 
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Integrate (2.5) over Q to obtain, 

dt [ p'^dx < {q-l)\\Vu\\Lo.^n) [ p'dx (2.6) 

i.e 

dtWphi < ^-^l|Vn||i.o(Q)||/)||L9 (2.7) 
Which imphes immediately 

\\p\\L4t)<C (2.8) 
with C independent of q, so our lemma follows. 

□ 

In the next proposition, we derive bound on norm of momentum pu. 
Proposition 2.2. Under condition (2.3), it holds that 

\\pu\\l^{Qt) ^ C:{\\po\\l^, \\uo\\lo^,\\V u\\li LOO, T), 0<T<T* (2.9) 

Proof. Let p be a fixed positive large number. Obviously, 

/ Po\uo\P+'^dx < c^Q^^ (2.10) 
Jn 

Without losing of generality, we assume p = 1. 

Multiplying \u\^u on both sides of the momentum equations in (1.1) yields 
1 d 



p + 2dt 



[ p\u\P+^dx+ I \Vu\'^\u\Pdx + p I \u\P{V\u\fdx 
Jn Jn Jn 



ap'^\u\^divudx + p / ap'^\u\^ ^u-\7\u\dx (^-H) 
n Jn 

= 1 + 11 

First, it follows from (2.4) and Holder inequality that 



|/| < C{ / \u\P\Vu\'^dxY''^{ / p\u\P+'^dx)^(^) 
Jn Jn 



V 



<- / \u\P\Vu\'dx + C{l p\u\P^''dx)^ 



2 



n Jn 



Similarly, 



II\<Cp{ \u\P\V\u\fdxy/'^{ p\u\P+'^dx)^(^ 



<? / l^r(V|n|)2 +Cp( / p\u\P+^dx)^ 



2 



n Jn 



(2.12) 



(2.13) 



A BLOW-UP CRITERION FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS 7 

Therefore, 



[ p\u\P+''dx + l [ \u\P\Vu\''dx + ?- [ \u\p(y\u\f 
^ ^ Jo. ^ Jo. ^ Jo. 



p + 2 

< C{p+l)i I p\u\P+^dx)^ 



dx 



n 



Integrating (2.14) over (0,T) yields 



/ p\u\P+\T)dx + - [ \u\P\Vu\'^dx + ^ I \u\P{V\u\fdx 

<Cip+l) fif p\u\P+^dx)l^^dt+^ [ polnor'dx 
Jo Jn P + 2 Jn 

< C{p+1){ I p\u\P+'^dxdt)^T^ + ^— I pQ\uQ\P+'^dx 
Jqt P + 2 Jn 



<C{p+l){[ p\u\P+'^dxdt)l^ + ^irc[;+2 



It follows that 



(2.14) 



(2.15) 



/ p\u\P-^\T)dx<C{p,CQ,T) 
Jn 

Moreover, it follows from Holder inequality and (2.15) that 

p\u\I^P+^Uxdt = [ {pl\u\P+^)l\u\P+^dxdt 
Qt Jqt 

< Cif pi\u\P+^dx)^ [ \u\^^P+^'^dx)Ut 
Jo Jn Jn 

<C Tif p\u\P+^dx)^[ |n|3(P+2))3dt 
Jo Jn Jn 

<Ci sup / p\u\P+^dx)i r {[ {\u\^fdx)-^dt (2.16) 
o<t<TJn Jo Jn 

<C{ sup / p\u\P+^dxf^ [ [ {Vlul'^fdxdt 
o<t<TJn Jo Jn 

<{C{p + 2f{[ p|^|P+2)^+^+2)| 
JQt 

<C{p + 2)f{[ p\u\P+'^dxdt)l +C{p + 2)f +Cl^'^^'^^ 

Jqt 



Here Co = 2co. Set 



5 10 

r = -,p + 2 = r^,C2 = —,ci = C{\pq\l^,\uo\l^,\\Vu\\lil^,T) (2.17) 
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We conclude from above the following reverse Holder inequality. 

[ p\uf^'dxdt<cir'^\[ p\ufdxdtY + cir''^'' + C^''^' (2.18) 

Jqt Jqt 



Define 



A{k)= p\uY''dxdt, /(A:) = cir^^fc^ 9{k) = (2.19) 
Jqt 

Then (2.18) could be written as 

Aik + 1) < f{k)A{kr + f{k) + g{kY (2.20) 

Write 

B{k) = max{A{k),g{k), 1), F{k) = 3f{k) (2.21) 
Without lose of generality, we assume that f{k) > 1. Hence 
B{k + 1) = maxiAik + l),g{k + 1), 1) 

< max{f{k)A{kY + f{k) + g{kY,g{kY, 1) 

< max{f{k)B{kY + f{k)B{kY + f{k)B{kY:g{kY, 1) 

< 3f{k)B{kY = F{k)B{kY 
By induction, we obtain 

B{k + l) < F{k)B{kY 

< F{k){F{k - l)B{k - lYY = F{k)F{k - lYB{k - lY'^ 
< 

< F{k)F{k - lYF{k - 2Y^ . . ■ F{2Y''~^B{2Y'' 

Consequently, 



B{k + l)rk+i < F{k)r'^'+^F{k-l)r'' ...F{2)^B{2)^ 



i=3 1=0 



(2.22) 



(2.23) 



(2.24) 



And 

fc+l k-2 

+ ^{k-i)r'-''-^ <oo (2.25) 



B{2) =max{A{2),g{2),l) =max{{ p\u f dxdt)T^ ,0^^%!) < oo (2.26) 

Jqt 

By definition 

( [ p\u\'''' dxdt)7^ < B{k)7^ < C (2.27) 
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Thus we can conclude , Vg' > 1, there exists a. k,a with < a < 1, r'^^'^ < q < 



Note that 



( / p'-'\ufdxdtp < ||p||^So'(Q^)( / p\ufdxdt)^ < C 
Consequently, 

\\pu\\l^{Qt) < C{\po\l^,\uo\lo^,\\Vu\\lil^,T) (2.28) 
Moreover, for a, /? > 0, p > 1 

im^hoo^Q^) < C{a, (3), im^hpiQr) < C{a, P,p) (2.29) 

□ 

The next lemma shows a connection between a convection term and the gradient of 
the density, which will play an important role in deriving the desired bounds on Vp. 

Lemma 2.3. Let F = put + pu ■ Vu. Then it holds that 

[ F^dxdt <C [ \Vp\'^dxdt + C, 0<T <T* 
Jqt JQt 

Proof. Note that 

/ F'^dxdt<C*{\\p\\Lo.(Q^)) pufdxdt + 2 \pu ■ Vu\'^dxdt (2.30) 

JQt JQt JQt 



ST WT 

It follows from proposition 2.1 that 



/ F^dxdt<C*{\\p\\Lo.(^Q^)) [ puldxdt + C{\\pu\\l^,Q .) f \yu\^dxdt (2.31) 
JQt JQt JQt 

Multiplying the momentum equation by ut and integrating show that 

j pufdx+ [ pu ■ Vu ■ utdx + -— j \Vu\^dx= j Pdivutdx (2.32) 



Note that 



/ Pdivutdx = dt Pdivudx - / Ptdivudx, (2.33) 
Jn Jn Jn 



and 

Pt + div(Pu) + (7 - l)Pdivu = 
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One gets 

/ Pdivutdx = 9t / Pdivudx+ / div(Pu)divudx + (7 - 1) / P(divu)^dx 

= dt Pdivudx- / (Pu) • Vdivudx+ (7 - 1) / P(divu)^dx 
Jo. Jn Jn 

Combing (2.32) and (2.34) yields 

— / |Vtipdx(r) + / pufdxdt + I pu ■ Vn • Ufdxdt 
2 Jn Jqt Jqt 

= \ I \yuo\^dx+ I Pdivudx(T)- / Podivuodx (2.35) 
2 Jn Jn Jn 

- Pu- Vdivudxdt + (7 - 1) / P(divu)^dxdt 

JQt JOt 



Direct estimates show that 



/ Pdivudx(T) <\l |Vupdx(T) + C (2.36) 
Jn 4 

/ pu -Vu ■ utdxdt < - I puf + 1 p\u ■ Vu\'^dxdt 
JQt 2 Jn Jn 



1 

< - 
- 2 



puj + C{\\p\u\'^\\loo(^q^)) I \Vufdxdt (2.37) 



T 



Q 



T 



puf + C 

On the other hand, it fohows from F = Au — VP and (2.2) that 

Pu ■ Vdivudxdt = / Pu ■ VdivA"^VPdxdt + / Pu • VdivA"^Fdxdt 



T 



. . (2.38) 

<C \Vp\^dxdt + e F^dxdt + C 
Jot Jqt 



Consequently, 

/ pu^dxdt + }- I \Vu\^dx{T) <C I \V p\^ dxdt + 2e [ F^dxdt + C (2.39) 
Jqt 2 Jn Jqt Jqt 



JT -Oil JQt 

Choosing e as 2C*e < 1, we may conclude 



/ F^dxdt<C [ \Vp\'^dxdt + C 

Jqt Jqt 

Which completes the proof of Lemma 2.3. 



□ 

We are now ready to show the desired L°°(0, T; L^(Q)) estimate of Vp. 
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Proposition 2.4. Under the assumption (2.1), it holds that 



sup 

o<t<r. 



f \Vp\^dx<C ,0<T<T* (2.40) 

/ pu^dxdt+ sup / iVnpdx < C ,0<r<r* (2.41) 
Jqt o<t<TJn 

|u||^2(f^)(it <C, 0<T <T* (2.42) 

Proof. Differentiating the mass equation in (1.1) with respect to Xj, 

dt{dip) + div(5ipu) + div(p9iu) = 

Which can be multiphed by 2dip to obtain 

dt\dip\'^ + div(|(9ip|2u) + |(9i/9pdivu + 29ipp9idivu + 2dipdiVi • Vp = (2.43) 

Integrating over 17 and using divu = divA^-^VP + divA^^F show that 

dt [ \d^p\^dx = - [ |a,ppdivudx-2 [ pdipdidivA~^VPdx 

- [ 2pdipdidivA-^Fdx- [ 2dipdiU-Vpdx (2-44) 
Jn Jn 

= -(^1+^2+^3 + ^4) 

Each term on the right hand side of (2.44) can be estimated as follows: 

\Ai{t)\ <\\divu\\L--it) / \dip\'^dx <\\divu\\Lo^{t) / \Vp\^dx (2.45) 
Jn Jn 

\A2{t)\ < C\\Vp\\l2\\VP\\l2 <C [ \Vp\^dx (2.46) 

Jn 

\A3{t)\<C\\Vph2\\Fh2 (2.47) 

|^4(t)| < C||Vn||Loo(t) [ \Vp\'^dx (2.48) 

Jn 

Consequently, 

dt [ \Vp\'^dx <C{\\Vu\\L--it) + l) [ \Vp\'^dx + C [ P'^dx (2.49) 
Jn Jn Jn 
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This, together with Gronwall's inequaUty yields 

/ |VH'd2;(t) <Ce^/o(llv«lk-(.)+i)ds( /■ |Vpo|2dx+ [\ [ F\s)dx)e-'' ^o(\\^-U^i-)+^)d- 
Jn Jn Jo Jn 

<C [ [ F^dxds + C 
Jo Jn 



<C [ [ \Vp\'^dxds + C 
Jo Jn 



(2.50) 



Hence 



sup 

0<t<T . 



I \Vp\^dx < C (2.51) 
Next, it follows from (2.31), (2.39) and (2.51) that 

/ pufdxdt+ sup / \Vu\'^dx<C (2.52) 
Jqt o<t<TJn 

This, together with Au = put + pu • Vu + VP, shows that 



<C + C\\Vu\\l^(^q,) + CWVpWl^q,) < c 



(2.53) 

□ 



Next, we proceed to improve the regularity class of p and u. To this end, we first 
derive some bounds on derivatives of u based on above estimates. 

Proposition 2.5. Under the condition (2.1), it holds that 

sup ||p^/2^i(i)|||2 + / \Vut\'^dxdt<C, 0<T<T* (2.54) 

0<t<T Jqt 

sup \\u\\h2 < C, < T < T* (2.55) 

0<t<T 

Proof. Differentiating the momentum equations in (1.1) with respect to time t yields 
putt + pu ■ Vut - Aut + Vpt = -pt{ut + u ■ Vu) - put ■ Vu (2.56) 

Taking the inner product of the above equation with ut in L^(r2) and integrating by 
parts, one gets 

d f 1 



dt 



f -pufdx+ [ \Vut\'^dx— j Pjdivutdx 
Jn 2 Jt^ Jn 

/ {pu ■ V[{ut + u ■ Vu)ut] + p{ut ■ Vu) ■ ut)dx 
Jn 



(2.57) 
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The last term on the left-hand side of (2.57) can be rewritten as (using (2.33)): 



Ptdivutdx = ^ / Jp(divu)2dx + / VP • (udivut)dx 
dt Jn 2 

+ J / {-Pu ■ V(divu)2 + (7 - l)P(divu)^)dx 
2 Jn 



(2.58) 



(2.59) 



It follows from (2.57) and (2.58) that 

Jt / (^P^? + |^(divu)2)dx+ / iVutpdx 

< / (2p|n| |nt| I VufI + p\u\ |nt| |Vn|^ + /9|n|^|nt| |V^n| + /9|u|^| | Vti| | Vuf| 
Jn 

+ p\ut\'^\Vu\ + \VP\\u\\Vut\ + 'yP\u\\Vu\\V'^u\ + j'^P\Vu\^)dx 

8 

Now, we estimate each Fi separately, where the Sobolev inequality and Holder in- 
equality will be frequently used. 

= / 2/>|ii||tij||Vnt|(ix 
Jn 

< C\\p'/ML-iQr)\\p'^'uth4'^uth2 (2-60) 

<e\\Vnt\\l.+C\\p'/\t\\l2 
Due to (2.40) and (2.41), one has 

sup / iViipdx < C 
o<t<T Jn 

Thus, it follows from Holder inequality, Sobolev imbedding and interpolation inequality 
that 



|-^2| = / /3|ti| I Vupfix 
Jn 

<c\\p^/ML^iQT) f p''^\ut\\yu\''dx 



< c\\p^/\t\\L4'^u\\l, 

^ ri\\ 1/2 ||1/2|| 1/2 ||1/2||V7 II IIV7 II 

< e||p^/V||L2||Vnt||i2 +C7||Vn||26 

< ellVutllia + C||pi/V|li2 + CllVull^a + CllV^nll^a 
<6||Vnt||i.+C7||pi/V|li2+Ch||^„ 



(2.61) 
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IF3I = / p\u\'\ut\\V'u\dx 
Jn 

<C||pVV||i.||V2^|U. (2-62) 
<C\\p^/^ut\\l,+C\\VMh, 

IF4I = / p\u\'^\Vu\\Vut\dx 
Jn 

<e\\Vut\\l,+C\\Vu\\l„ 

l^sl = / p\ut\^\Vu\dx 
Jn 

<C\\puUl4'^u\\l2 



< C\\p'/'u,\\% 



in 

<C\\Vu\\l2\\vMl^, 
Finally, noting that Vu G L'^{0,T; L'^{n)) n L^{0,T; L^{n)), 

L^{Qt) C L°°(0, T; L'^in)) n L\0, T; L'^in)) 



Hence 



\Fs\ = [ -f^P\Vufdx 
Jn 

<C [ \Vu\^dx 
Jn 

<C\\Vu\\l^^u) [ |V 
Jn 

< C||Vu||ioo(Q), 



u\ dx 



(2.63) 



(2.64) 



<e||pV\,||i«+C||pVV||i. 
<e\\ut\\l,+C\\p'/^u4l, 
<e\\Vut\\l, + C\\p'l'uA\l,, 

iFel = / \VP\\u\\Vut\dx 
Jn 

<C [ \Vp\\Vut\dx (2-65) 
Jn 

<e\\Vnt\\l, + C\\Vp\\l„ 
IF7I = / -fP\u\\Vu\\V^u\dx 

(2.66) 



(2.67) 
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Collecting all the estimates for Fj, we conclude 



(2.68) 



Jt jjh^''^' + |p(divu)2)dx + ^ iVutl^dx 

< 5e / \Vut\'^dx + C(||p^/V|||2 + \\u\\]j2 + \\Vp\\l2 + \\Vu\\l^) 
Jn 

Thanks to the compatibility condition: 

11 1 1 

poix)2{pQ{x)2ut{t = 0,x) + p§uo- Vuo{x) - p^g) = (2.69) 

it holds that 

Po{x)^ut{t = 0,x)= pIuo • Vuo{x) - pig G L\Q) (2.70) 
Therefore, for arbitrary small e, (2.68) yields 

sup ||p^/2^t(t)|||2 + / \Vut\'^dxdt<C, 0<T<T* (2.71) 

0<t<T Jqj. 

Moreover, 

Au = put + pu-Vu + VP € L'^L^ 



Hence, 

""■'1^2 



0<T<T* 

Our lemma follows immediately. 



sup ll-ull^a < C (2.72) 



□ 



Finally, the following lemma gives bounds of derivatives of the density and the second 
derivatives of the velocity. 

Lemma 2.6. Under the condition (2.1), it holds that 

sup (llpt(t)llL^o + WpWw^m) <C, 0<T<T* 
0<t<T 

fT 

/ Mt)\\l,2.,odt<C, 0<T <T*,qo = min{6,q) 
Jo 

Proof. It follows from (2.71) and (2.72) that 

ut G L\0,T;L^{n)),Vu G L%Qt) 

F G L^{0,T;L^{n)) 
Differentiating the mass equation in (1.1) with respect to Xi, one gets 

dt{dip) + div(5i/9u) + div(p9iu) = 
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Multiplying above identity by qQ\dip\''°~'^dip, one gets that for some go = min{<6,q), 
dt\dip\'i'' + d\Y{\d,prn) + (qo - l)|aip|'i°divu 

(2.73) 

+ qo\d^p\'i'^-^dippd^diwn + qo|ai/)|^o-2g.^^.^^ • Vp = 
This, together with divu = divA^^VP + divA~-'^F shows that 

dt [ \dip\^°dx = -{qo-1) [ l^ipl'^odivudx-qo / pldipl'^^'-^dipdidivA-^VPdx 
Jn Jn Jn 

- [ qop\^ip\'^'>-^^^p^^divA-^Fdx- qo [ \dipr-^dipdin ■ V pdx 

= -(^1 + 52 + ^3 + ^4) 



Which can be estimated as: 



\B3{t)\<C\\\Vpr-^\\ ^||F||l,o 
<C\\Vp\\f-'\\FU.o 



(2.74) 



\Biit)\<{qo-l)\\Vu\\L--it) [ \dip\'^°dx<C\\Vu\\L--{t) [ \Vp\'^°dx (2.75) 

Jn Jn 

|S2(t)| < C|||Vp|*^-^|| ^||VP||l,o < C / |V/)|«°dx (2.76) 
iw-i J^ 



(2.77) 



\B4it)\<C\\Vu\\Loo{t) [ |Vp|9odx (2.78) 

Jn 

It follows from (2.74) - (2.78) that 

dt\\Vp\\L'>o < C(||Vtx|Uoo(t) + l)||Vp||L,o +C||F||l,o (2.79) 

Hence, 

sup ||V/?||£90 < C 

0<t<T 

Therefore, due to this, (2.71) and interpolation inequality, one has 

p^ = -(^u■Vp + pdivu) G L'^L^" (2.80) 

Recall that 

Au = F + VP £ L^L'i^ 
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gives 



T 







\u\\l,2„o^^)dt<C (2.81) 



This will close the estimates and guarantee to have an extension of the strong solution. 

In fact, in view of Proposition 2.4 — 2.5 and lemma 2.6, the functions {p,u)\t=T* = 
limt^T'ip, u) satisfy the conditions imposed on the initial data (1.7) — (1.8) at the time 
t = T* Furthermore, 

p^ut + p^u ■ Vu G L°°L^ 

— Au + VP\t=T*= lim P^ip^ut + P^u -Vu) = p^g\t=T* (2.82) 

Where g\t=T-' G L^i^). Therefore, we can take (/?, u)\t=T* as the initial data and apply 
the local existence theorem [T3] to extend our local strong solution beyond T*. This 
contradicts the assumption on T*. 

□ 

Note that after some minor modifications, the above ideas also works in both periodic 
case and = i?^, so theorem 1.2 holds. 

When Q = T^, using an estimate from Desiardin[10j. 
rT 



I 



p2Ut\\^2('r2) + WP^U ■'^U\\^2(rp2))ds + sup ||V'u(t, )||^2(r2) 

o<t<T ^ ' 



< exp(Cexp( / \\p\\iao/j'2\5{s)ds)) 
Jo 

Where 

6{s) = 1 + \\Vu\\l2^t2) G L^{0,T) (2.84) 
it follows from (2.83) - (2.84) that 

/ F'^dxdt+ sup \\Vu\\l2(rp2)<C (2.85) 

jQrp 0<t<T ^ ' 

Then we can do a similar estimate step by step as three dimensional case to obtain a 
higher regularity of {p, u) . We omit the detail for simplicity. 
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